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Abstract 

Using a recent next-to- leading-order calculation of the photoproduction dou- 
ble differential cross section for heavy quarks, we study the possibility of ex- 
tracting the gluon density of the proton from heavy-quark photoproduction 
data. We discuss the theoretical uncertainties connected with this method, 
and we conclude that they are well under control in a wide x domain. 
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The gluon density of the proton is usually extracted from fits to deep inelastic 
scattering data. There it affects the evolution of the quark density, according to the 
Altarelli-ParisiB evolution equation. Its determination is indirect, and to some ex- 
tent, it is influenced by the particular parametrization chosen for the fits. A direct 
determination is therefore highly desirable. In hadron-hadron collisions, processes 
such as direct photonffl and jet® production are directly sensitive to the gluon dis- 
tribution. Although these processes provide further constraints on the gluon density, 
they do not allow a precision comparable with the DIS experiments. 

The ep collider HERA offers new opportunities to measure the gluon density 
directly. Since the initial state is simpler than in hadron-hadron collisions, the gluon 
density enters in a simpler way. Furthermore, various methods can be used, thus 
allowing for consistency checks. Some of the methods are reviewed in ref. ||. The 
most promising one is based upon the measurement of the longitudinal structure 
function® . It will probably allow a determination of the gluon density in the x range 
between 1(T 3 and 5 x 1CT 2 for 25 GeV 2 < Q 2 < 150 GeV 2 . 

There are a number of methods that rely upon the study of hard production pro- 
cesses. At HERA, these processes arise both from the pure photon component and 
from the hadronic (sometimes called resolved) component of the photon. In jet pro- 
duction, the effect of the latter is strongly enhanced, because the typical elementary 
process gg — > gg is about 20 times larger than the processes jg — > qq and -yq — ► qg. 
One possible strategy to suppress the hadronic component is to consider only off-shell 
photon jet production. This suppresses the rate, but eliminates the hadronic com- 
ponent. One should then disentangle the two pure photon processes jg — > qq and 
75 — > qg. The other possibility is to require heavy-quark production. In the hadronic 
component only the two subprocesses (at leading order) gg — > QQ and qq — > QQ 
contribute, and they are two orders of magnitude smaller!! than the typical process 
99 99- The hadronic component is not eliminated completely with this method, 
but there are large kinematical ranges in which it is negligible. 

Along this line a popular method is to consider J ftp production H . Although 
experimentally very convenient, this method suffers from several theoretical ambigu- 
ities, essentially due to difficulties in the computation of the J/ip production cross 
section. Open-charm production is a more promising method®. It is less affected by 
theoretical ambiguities, and it is sensitive to the gluon density in the small- x region. 
In ref. |§ a study of the feasibility of measuring the gluon density using open charm 
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production was performed, using leading-order cross sections. It was concluded that 
f&)(x) can be measured for x in the range between 10 _1 and 10 -3 . 

In order to extract the gluon density from open-charm production to next-to- 
leading accuracy, a next-to-leading calculation of the heavy-quark photoproduction 
cross section is needed. Such calculations for single-inclusive cross sections (i.e. inte- 
grated over the momentum of one of the heavy quarks) have been available for some 
time0'0. They are not, however, very useful for our purpose, since they do not 
allow for the full reconstruction of the x variable of the parton in the proton from the 
final-state observables. 

We have recently completed a next-to-leading calculation of the double differential 
heavy-quark photoproduction cross section. This calculation, which will be described 
in more detail in a forthcoming publication^! , is an extension of the results presented 
in refs. [[y|] and [14]. In the present work, we will apply this result to the problem of 
extracting the gluon density of the proton from open-charm production. The main 
goal of this work is to show that the theoretical uncertainties associated with this 
procedure are well under control, and to assess their magnitude. The possible sources 
of uncertainty in the problem are the following: radiative corrections of even higher 
order (i.e. C(a em a|)), not included in our calculation, the poor knowledge of the 
value of the heavy-quark mass, uncertainties on the value of Aqcd, and the presence 
of a hadronic component of the photon. Furthermore, one should not forget the large 
hadronization effects that are usually found in charm production. There are good 
reasons to believe that these effects are smaller in the photoproduction than in the 
hadroproduction case. We will also argue that there are several ways of getting rid 
of the major hadronization ambiguities. 

We begin by writing the heavy-quark cross section at the leading order in the 
following form 



dy^dM 2 9 dx g dM 2 Q73 E 2 



where M Q q is the invariant mass of the heavy-quark pair, and y Q Q is the rapidity 
of the pair in the electron-proton centre-of-mass frame (we choose positive y in the 
direction of the incoming photon). E is the electron-proton centre-of-mass energy, 
and 

x 7 = ^pexp(y Q p) (2) 
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x a = ^^ ex P(-2/Qg)- (3) 

The function ffl is the photon density in the electron in the Weizsacker- Williams 
approximation, 

A , )(x ^ l) = ^l + i^l log ^0l (4) 

and f^' is the density of gluons in the proton. The Born level partonic cross section 



(tQ){s) is given by 



e 2 Q a em a s (p R ) nf3p 



3-f3\ l + p 



log - -4 + 2/3 



(3 °l-(3 



19 V / m 2 4 

where m is the heavy-quark mass, e Q is its charge in electron charge units, and 

4m 2 



(5) 



P 



(3 = s/l^~p . (6) 



s 

The renormalization scale is set to fi R = /x , and the reference scale p is defined by 



//o = sj(f T +f T )/2 + m\ (7) 

where j»t and pT are the transverse momenta of the heavy quark and of the heavy 
antiquark. The factorization scale for the proton is set to 2p . 

Assuming that we identify the left-hand side of eq. ([I]) with the data, we can invert 
the equation, and get a first determination of fjP': 

J^data rp2 

ho) ( U 2\ = X _EZ £: (8 ) 

Jg \ gi r'Fl x 9 j j M 2 .(e), 2 \ - (0) i *jT2 \' K ' 

The inclusion of radiative corrections does not pose any problem. We write the full 
cross section as 

X9 dx g d d^ = h #^^)#^ (9) 

where A represents all the radiative effects. In A we have also indicated explicitly 
the functional dependence upon the gluon density in the proton. The light quarks, 
which enter at the next-to-leading order via the ^yq — > qQQ process, give a small 
contribution (less than 5% for all values of x g and M Q q considered here). The A term 
is given by 

A(fW,x g ,M 2 Q7 ;) = ^fdxdx^^x^^f^ix^^a^aUpn) 

3 

da?} 



x ■ 



dM^dy Q Q 



X^xE , Mqq, ijQQi Pfi Pri P^)i (10) 
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where 

1 M 2 - 

jLi tXj t\j g J— J 

is the heavy-quark pair rapidity in the partonic centre of mass frame. The factoriza- 
tion scale of the photon /x 7 is set to 1 GeV (we verified that the results are rather 
insensitive to the choice of /i 7 ). For a complete discussion of all the scale choices we 
refer the reader to a forthcoming publication^ . 

We now write the full as 

#W) = /f fcAO + tfW)^ (12) 

where the second term is the next-to- leading correction, and plug it back into eq. 
We get 

E*A(f(°\x 9 ,M? 
#(x 7)M g)4?(M^) 



Jo { X 9^f) - „( e ), 2\ -(0)/— ' ^' ^ ' 



We have neglected the piece contained in the A term, the corresponding contri- 
bution being of order a em a|. It could also be easily incorporated by iterating eq. (|T3D , 
using the full structure function in the right-hand side. In the simple illustration we 
have given, we have integrated the cross section over all variables but x g and M Q q. 
In general, in realistic experimental configurations, other cuts will be applied to the 
data. The procedure for the extraction of outlined above can be carried out also 
in this case, provided the same cuts as applied to the hadronic final states in the 
data are also applied to the partonic final state in the calculation. One interesting 
possibility is to apply this procedure to the invariant mass of the two-jet system, 
instead of the invariant mass of the heavy-quark pair. This makes no difference at 
the Born level, but next-to-leading corrections will turn out to be different, and one 
should compute them with an appropriate jet definition (which should match the jet 
definition used in the analysis of the data) in order to get a meaningful answer (see 
refs. [jy| for a discussion of this point). Furthermore, with this procedure there would 
be no need to identify both charmed mesons, and the result would be less sensitive 
to fragmentation effects. 

We now begin with a study of the charm differential distribution in the variable x g 
(as defined in eq. (j|)) at HERA, for various cuts in M Q q. These are the distributions 
which are relevant to the extraction of the gluon density at HERA. All our results 
are obtained for v^S* =314 GeV. In fig. [I] we plot the leading-order and next-to- 
leading-order pure photon cross section (throughout this work, we use the proton 
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structure function set HMRS B of ref. |l(|). From the figure we first notice that 
radiative corrections are moderate, but not negligible. This is an indication that 
the perturbative expansion is reliable in this case. We also observe that the effect 
of the radiative corrections cannot be described by a simple if-factor. From fig. [TJ 
neglecting for the moment the contribution of the hadronic component of the photon, 
we can infer that the distribution is sensitive to the gluon density of the proton down 
to x g of the order of a few times 10~ 4 . 

We now examine the effect of the hadronic component of the photon. This con- 
tributes another term to the cross section, given by 



da H = J2[ dxdx, If' f^(z,^) (2>,iA - 

■ ■J Jx 1 \ Z / Z 



fPfafj&d&ij. (14) 



This is formally of the same order as eq. ([]]), since the photon parton densities 
are of order a em /as- The contribution of the hadronic component, including next- 
to-leading orders (computed using the program of ref. Hl3| ) is displayed in fig. || for 
two different parametrizations of the photon structure function. We have chosen 
the set LAC 1 (ref. |7|) and the set ACFGP HO-mc (ref. @). These two sets 
have extreme gluon distributions (as can be easily seen from fig. 10 of ref. |19fl ), 
and they may therefore give an idea of the magnitude of the hadronic component 
contribution (experimental results from HERA will help to better specify the value 
of this component). We see from fig. |2| that the hadronic component gives a large 
contribution for the smallest invariant-mass cut. Even in this case, however, there is 
a region of small x g in which the hadronic component is negligible with respect to 
the pure photon term. When we increase the invariant mass cut, we notice that the 
hadronic component decreases faster than the pure photon contribution. This is due 
to the fact that in the hadronic component, when we increase the invariant mass, the 
production process is suppressed by the small value of the gluon density of the photon 
at large x. We therefore see that a large region in x g is accessible by this method. 
By pushing the invariant-mass cut to 20 GeV we can reach the region of x g = 10 _1 . 
Observe that statistics will not be a problem even at these large invariant-mass cuts. 
With 100 pb _1 of integrated luminosity, there would be about 10 5 events per bin in the 
case of a 20 GeV cut (before accounting for experimental efficiencies). We therefore 
find that the conclusions of ref. || are not spoiled by theoretical ambiguities due to 
higher-order effects or to the hadronic component of the photon. We then conclude 
that the theoretical uncertainties in the heavy-flavour cross section, in the range of 
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10 3 < £ 9 < 10 l , are small enough to allow for a determination of the gluon density 
in the proton. 

The authors of ref. reach a conclusion that contrasts with ours, that is to 
say that the ambiguities related to the hadronic component of the photon spoil the 
predictivity of the method. This conclusion is based upon a single inclusive calcula- 
tion. They find that in particular the set LAC 3 of the photon structure functions 
of ref. ]17| gives a hadronic contribution that competes with the pure photon contri- 
bution for all rapidities. We find that even with the LAC 3 set, when looking at the 
double differential quantities we have considered, it is still possible to perform the 
extraction of the gluon density in the proton, although in a more restricted x g range. 
However, we have chosen not to include this structure function set in our analysis for 
the following reasons. The LAC 3 set has an unphysically hard gluon structure func- 
tion. For Q 2 = 5 GeV 2 , it peaks at x = 0.9 and carries 70% of the total momentum 
of the hadronic component of the photon. Because of this very pronounced peak, it is 
roughly as hard as the pure photon. We have examined the origin of this behaviour. 
As can be seen from ref. ]T7|], the set LAC 3, unlike the sets LAC 1 and LAC 2, 
is obtained by fitting data for F^x^Q 2 ) at Q 2 values down to 1.31 GeV 2 (for the 
other sets the fits start at Q 2 = 4.3 GeV 2 ). From Q 2 = 1.3 GeV 2 to Q 2 = 4.3 GeV 2 , 
F%{x, Q 2 ) grows rapidly in the region x > 0.2. In order to reproduce this growth us- 
ing the QCD evolution, one is forced to assume a very hard gluon component, which 
under evolution will feed down a large quark component at moderate values of x. On 
the other hand, at this low Q 2 values, there may be other reasons, having nothing to 
do with the QCD evolution, that can cause the observed growth. For example, new 
thresholds for single resonances or resonance-pair production can be opened. After 
all, at Q 2 = 1.3 GeV 2 and x = 0.3 the mass of the produced hadronic system is only 
of 1.74 GeV. We therefore believe that these data are more consistent with a change 
in regime for from a low-energy VMD behaviour to a high-energy point-like be- 
haviour, than with perturbative evolution; one should thus not attempt to fit it with 
QCD evolution alone. Indications that one should use Q 2 values larger than 4 GeV 2 
were also given in ref. Furthermore, the authors of ref. [17] also express some 

scepticism with respect to their set 3. Finally, experimental results on jet production 
in photon-photon collisions!! and photon-hadron collisions! 2 ! clearly disfavour the 
set LAC 3. 

We will now turn to a discussion of the other theoretical uncertainties involved 
in our procedure, in order to estimate the precision of the method. One source of 
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uncertainty is due to the truncation of the perturbative expansion. We estimate this 
uncertainty by varying the renormalization and factorization scales by a factor of two 
above and below their reference value. The result is plotted in fig. |3|. We see that the 
effect of the variation of the factorization scale is moderate, while the renormalization 
scale uncertainty amounts to an uncertainty of ±10 to 20% on the result. 

In fig. [| we show the dependence of our result upon the heavy quark mass. Even 
in this case the variation of the cross section is of the order of 10 to 20%. We should 
remind, however, that the mass of the heavy quark (unlike the renormalization scale) 
is a physical parameter that enters the cross section. It therefore makes sense to reduce 
this uncertainty by using values of the mass that give a good fit of the data. There 
are of course also uncertainties associated with the error on the present determination 
of Kqcd- In practice, it will be more convenient to measure asfg(x), which is much 
less sensitive to the value of Kqcd- 

We observe that perturbation theory alone makes a prediction for the ratio of our 
differential distributions for different mass cuts in the x g region in which the hadronic 
component is small. In leading order and up to scaling violation, the gluon density 
in the proton cancels out in this ratio. We present a plot of these ratios in fig. |5]. As 
can be seen from the figure, the mass and scale dependence of the ratios for the pure 
photon contribution is negligible, while significant changes can be observed when the 
hadronic component is included. This in principle might allow for an independent 
check of the hadronic structure of the photon. 



As suggested in ref . , an additional help to separate the pure from the hadronic 
photon contributions comes from tagging the photon energy by measuring the energy 
of the recoiling electron in the small-angle luminosity monitors. This approach might 
have, in general, problems at next-to-leading order: processes such as — > qQQ, 
the light quark being emitted preferentially in the direction of the photon, will look 
and will be reconstructed as hadronic photon events. As was mentioned previously, 
however, the contribution of these processes is numerically negligible. The photon- 
tagging technique could then help to further constrain the gluon density of the photon. 

We conclude by giving in table 1 the total cross sections for b and c production at 
HERA. Observe that the sensitivity to the variation of the scale and mass parameters 
is much smaller for b than for c production. This suggests that by using b production 
a better precision may be achieved in the determination of the gluon density. On the 
other hand the x g region covered is smaller. We see in fact that using the structure 
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function set MSRD- (a set with a gluon density more singular at small x) the charm 
cross section grows by 80%, while the b cross section grows by less than 10% with 
respect to our central value. 

Observe the large difference in the hadronic component when using the two differ- 
ent sets of photon structure functions. When going from the ACFGP to the LAC 1 
set, the hadronic component of the charm cross section grows by a factor of 9, while 
for bottom it grows only by a factor of 3. This difference is mainly due to the small 
x growth of the gluon density in the photon. It is likely that the charm cross section 
and distributions at HERA will also help to constrain the gluon density in the photon 
at small x. 
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a C c (nb) 


a b i (nb) 


Central value 


423 


4.37 


fi R = 2/i 


362 


3.90 


H R = /jlq/2 


489 


4.94 


Hf = 2^o 




4.50 


A*f = Aio/2 




4.13 


Low mass 


705 


5.16 


High mass 


270 


3.72 


A 4 = 100 MeV 


337 


3.66 


A 4 = 300 MeV 


453 


4.66 


MRSD- 


765 


4.65 


Hadr. comp. LAC 1 


723 


3.10 


Hadr. comp. ACFGP 


80.5 


0.91 



Table 1: Charm and bottom pair-production cross sections at HERA. The pure 
photon component and the hadronic component are shown separately. The central 
values for the scales are fi R = m c for charm and m& for bottom, /i F = 2m c for charm 
and rrib for bottom. The central mass value for charm used here is 1.5 GeV, the low 
value is 1.2 GeV and the high value is 1.8 GeV. For bottom we set mb = 4.75 GeV, 
the low value is 4.5 GeV and the high value is 5 GeV. The structure function set for 
the central value is HMRSB, with A 4 = 190 MeV. When studying the sensitivity to 
A we use two HMRSB fits@ obtained with the values of A indicated in the table. 
We also show the values obtained using the MRSD- structure functions, which have 
stronger enhancement at small x. For the hadronic component only the central values 
are used, since the uncertainty due to the choice of the photon structure function is 
by far the largest. 
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Figure Captions 

Fig. 1: Differential cross sections for charm production, histogrammed in the log- 
arithm of x g = exp(—y Q Q)M Q Q/E, for m = 1.5 GeV, HMRS B structure 
functions for the proton, and default values for the scale choices. Only the 
pure photon component is included. 

Fig. 2: Contribution of the hadronic component of the photon, compared with the 
pure photon component, for two extreme choices of the photon structure 
functions. 

Fig. 3: Scale dependence of the cross section (pure photon only). The solid line 
corresponds to the choice ji R = fi and \i F = 2/i , the dotted line to \i R = 
yUo/2 and fi F = 2/i , the dashed line to fi R = 2/i and ji F = 2/i , and the 
dot-dashed line to jj, R = fj, and fi F = /U /2. This last curve is not shown for 
the smallest mass cut, because the scale goes outside the range of validity 
of the structure function parametrization. 
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Fig. 4: Charm mass dependence of the cross section (pure photon only). The 
solid line corresponds to the choice m c = 1.5 GeV, the dotted line to 
m c = 1.8 GeV and the dashed line to m c — 1.2 GeV. 

Fig. 5: Ratios of cross sections for different invariant-mass cuts, (a) Compari- 
son between the pure photon case and the full (pure + hadronic) one, for 
two different photon parton distributions, (b) Same as (a), for a different 
invariant-mass cut. (c) Mass dependence of the ratio in the pure photon 
case, patterns as in fig. 4. (d) Scale dependence of the ratio in the pure 
photon case, patterns as in fig. 3. 



